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C~^- ' In supersymmetry with large tan j3 the decays B°(B°) — ► l + l arc dom- 

inated by the scalar and pseudoscalar Higgs penguin diagrams leading to 
strong enhancement of leptonic decay rates with potentially large CP asym- 
metries in the r + r~ decay modes measurable in BELLE or BABAR exper- 
. iments. The TAUOLA r-lepton decay library supplemented by its universal 

interface can efficiently be used to search for B°(B°) — > r + r~ decays, and 
to investigate how the CP asymmetry is reflected in realistic experimental 
observables. 
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1. Introduction 

Understanding the origin of CP violation is one of the most important 
tasks of particle physics. Until now, CP violation has been firmly established 
in K- and -B-physics in a series of high statistics experiments. One way 
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of testing the conventional CKM description of CP violation is by more 
precise measurements and joint analysis of CP asymmetries measured in 
different channels with a hope of finding some inconsistency signaling a 
contribution from physics beyond the Standard Model (SM) to CP violation. 
Alternatively, observation of a non-zero effect in channels in which no (or 
negligibly small) effects violating CP are predicted by the SM would be an 
unambiguous signal of new physics contribution to CP violation. 

In this talk presented are results of recent analyses of flavour changing 
decays of the neutral B mesons into lepton pairs, B® — > l + l~ m, for which 
the SM predicts no CP violating effects. This decay mode is very sensitive 
to new physics which affects the 6-quark Yukawa couplings j^l El |I] . Ap- 
proximate and full one-loop calculations in the supersymmetric extension of 

* Talk given at the International Conference "Matter to the Deepest", Ustrori, Poland, 

8-14 September 2005 
t Speaker. E-mail address: Jan.Kalinowski@fuw.edu.pl. 
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Mode 


SM expectation 


Exp. limit Ref. 


B d -» |U+// - 

B d -> r+r- 
S s -> r+r- 


1 x l(T iU 

2.8 x 1(T 8 
3.7 x 1(T 9 
1 x 1(T 6 


8.3 x 10~ s [S] 
3.8 x lfT 8 [| 
2.7 x 10" 3 P3J| 
1.5 x 1(T 7 H 
~ 5 x 10 -2 [HI 



Table 1. Upper limits for the branching fractions of leptonic B- meson decays. 

the SM with large tan (3 (the ratio of the vacuum expectation values of the 
two Higgs doublets) EJ E] showed that truly spectacular enhancement of 
the rates of the decays — » can be expected. Moreover, new physics 
can also lead to observable CP violation in these decays pQ. CP violation 
could manifest itself through non-equal leptonic decay rates of the B°(t) and 
B°(t) states (tagged t = as B° and B°, respectively). If polarization of 
final state leptons can be determined, additional information on CP viola- 
tion could be provided by non-equal T(B°(t) — > and T(B°(t) — > IrIr) 
[or T(B°(t) -» Z^Z^) and T(i? (t) -» Ipl)] decay rates HE]. 

Leptonic decays are theoretically clean as the only non-perturbative 
quantities they depend on, are the B° meson decay constants Fs d s , which 
cancel out in suitably defined CP asymmetries. However, none of these de- 
cays have been seen so far: the best upper limits on the B^ s — > and 
t + t~ branching fractions at present are listed in Tabled 

With the rates as predicted by the SM, the detection of the n + fJ>~~ de- 
cay channel will become possible only at the LHC; in the hadronic collider 
the t + t~ channel is extremely challenging. New physics (like supersym- 
metry) can increase significantly their rates to a level that they can be 
observed at BaBar, BELLE or Tevatron in near future. We find, however, 
that the ratio of time integrated leptonic decay rates of B°(t) — ► fJ> + fi~ and 
B°(t) — > is unlikely to deviate appreciably from unity. Polarization 

measurement also seems very difficult in the case of the /i 4 ~ \x~ channel. In 
the t + t~ channel the situation can be quite different: large CP violating 
effects can be expected, and r polarization measurement is possible. Re- 
sults of first experimental analyses of this channel have recently been made 
public jlUj . Increasing the experimental efficiency of identifying r leptons 
and suppressing the background might result in a much stronger limit or, 
hopefully, in observing the signal in B meson factories. We show that the 
existing TAUOLA package |12j and its universal interface |13j may prove 
useful to search for these decays and for the CP violation since in realistic 
scenarios of large tan (3 MSSM the decay rates can significantly be increased 
to a level measurable at the running BaBar and BELLE experiments and 
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the CP asymmetry in the B®(B®) — > t + t~ channel can be quite large and 
potentially measurable. We also identify, in addition to the ratio of time 
integrated leptonic B d (t) and B®(t), two realistic CP-sensitive experimen- 
tal observables which reflect r-lepton polarisations: the ir^ energies from 
r — > irv decays and the acoplanarity angle between the decay planes of the 
p mesons which originate from r — > pv. The former is sensitive to the longi- 
tudinal, while the latter to the transverse polarizations of r's coming from 
B® and B® decays. 

2. Preliminaries 

In leptonic i3-meson decays if, for example, 

T(B°^lil L )^T(B°^lp R ) , 
T(B° - l+lx) ? T(B° - f+li) (1) 

CP is violated because the initial and final states on both sides transform 
into each other under CP |Sj • The amplitudes of B° decays into two helicity 
eigenstates read 

A L = (ltl L \B ) =M B (a + b (3) , 

A R =(lp R \B°) = M B (a-bp) , (2) 

Similar formulae with a and b replaced by a and b, respectively, give the 
amplitudes Al and Ar for the corresponding B° decays. Here (3 = (1 — 
Amf/Ml) 1 / 2 and o, b, a and b are the coefficients in the effective Lagrangian 
describing B°(B°) — ► l + l~ decays 

£ eff = Bl d ^(b Std + a s>d j 5 )ipi + B°Jji(Ps4 + a-.dT 5 )^ • ( 3 ) 

(the subscripts d and s referring to non-strange and strange B° mesons, 
unless explicitly written, will be omitted). Hermiticity (CPT invariance) 
implies b = b* and a = —a*. 

Since in leptonic decays no strong phases are involved, inequality (^Q) can 
occur only through the mixing of the B° and B° mesons. In the standard 
formalism |14| the state which at t = is a pure B° (B°) evolves in time as 

\B° phys (t)) = g+(t)\B°) + - p g-{t)\B Q ) (4) 

(for Bph ys (t) replace B° <-» B° and p q). The coefficients g±(t), neglect- 
ing the difference of the decay widths of the two B° mass eigenstates and 
denoting AM = M B o - M B o <C M B = (M#> + Af B o)/2, read 

g+(t) = e - lMBt -^ cos — t , <?_(*) = e -^'-2* i S in — i , (5) 
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and the ratio p/q is calculated from the effective Hamiltonian 

p/q = (HyH 12 f/ 2 (6) 

with H 12 = Mi 2 + jr 12 = {B°\H cS \B°), etc. M- 

Since \A L \ = \Ar\, \Ar\ = \A L \ and A L /A L = (Ar/Ar)* , CP is vio- 
lated if either 

\q/p\ ? 1 or Im(A L ) ± Im(A^) , 

where Xl = qAi/pAL and Xr = qARjpAR . 

The simplest quantitative measures of CP violation are provided by the 
asymmetries constructed out of time integrated polarized decay rates 



fJl dt [r(B° phys (t) -> Ijlj) - r(ff p ° hys (t) -> l + R l- R )] 



u III dt [r(5 p ° hys (t) - i + R r R ) + r(B° phys (t) - ip L )] • 

and the ratio of integrated unpolarized decay rates 

J? * r(B° hys (t) - /+/-) 



If the statistics of tagged events is low, or experimental determination of the 
decay time t is difficult, the fully integrated asymmetries Aq P = Aq P (0, oo), 
Aqp = Aq P (0, oo) and the ratio Ri = Ri(0, oo) can be exploited. If \q/p\ = 1 
(as in the SM and many SUSY extensions), the expressions for the asym- 
metries Aqp, Aqp and the ratio Ri simplify to 

! -2 x Im Xl . 2 -2 a ImA R 

^CP " 2 + x 2 +x 2 | Al |2 > ^CP " 2 + x 2 + x 2 \ Xr \2 ' ^) 



(|^l| 2 + |^| 2 ) (1 + x 2 ) - x { |^| 2 Im(A L ) + i^j 2 Im(A fl )} 
(|^l| 2 + W) (1 + x 2 ) + x{\A L \nm(X L ) + |^| 2 Im(A fi )} ' 



where x = AM/r. The asymmetries Aq P , Aq P , as functions of Xl,r, are 
therefore bounded from above by [H] 



A 1 ' 2 
^CP 



<{2 + x 2 )- 1 ' 2 . (12) 



which has important consequences. Since x s > 20.6 for the B®-B® system, 
the CP asymmetries in the leptonic B®(B®) decays can reach at best ~ 4.5% 
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irrespectively of the amount of CP violation. In contrast, for the B®-B® 
system, for which Xd = 0.771 ± 0.012 they can be as large as ~ 60% and 
hopefully measurable at BaBar and BELLE in a relatively clean environ- 
ment. Therefore below we will consider only the CP asymmetries in the 
B° d {B° d ) -» l+r decays. 

The time dependent i3°-meson mixing can easily be dealt with by in- 
troducing time dependent effective factors a e s , b e s , a e ff and 6 e ff • The factor 
a c fj is defined as 

cieff(t) = a g + (t) + a - g-(i) (13) 
V 

and the other factors are defined in a similar manner jT]. Then the instan- 
taneous B° widths into left- or right-handed I's read 

T(B° phys (t) -» l+ll [1+1$) = ^-0 loeff(t) + [-]/3 6 eff (t)| 2 (14) 

and those for B° are given by similar formulae with a e s(t), b e g(t) replaced 
by a e fr(t), b e fi(t). CP is violated because in general a e fi(t) / — and 
b eS (t)^b* eS (t). 

In the case of the t+t decay mode the r polarisations are best identified 
by measuring the tt+ energy spectra from r — ► ttu decays |15j . The spin 
density matrix formalism of ref. ^Hl allows to construct also observables 
sensitive to transverse polarization of the final state r's. The spin weight 
for the complete event (B — > tt —* decay products) can be written in the 
form 

WT=^(1+ £ R H h \ 4 + E R * h\ + £ R 0J 4) • (15) 

ij=x,y,z i=x,y,z j=x,y,z 

where the polarimetric vectors hi and /12 are determined solely by the dy- 
namic of the r decay processes, and 



Roo — 


+ 1, 


1 I — Ryo — Rox 


— Roy 


= 0, 


Rzz = 


-1, 


Rxz — Ryz — Rzx 


— Rzy 


= 0, 


Roz = 


—Rzo = 


2Re(a cfr 6* ff ) 

l^efFl 2 /? 2 + |«efr| 2 ' 






Rxy = 


Ryx 


2Im(a eff & e * ff ) 


Rxx 


B 

— ^yy — 


|6 eff | 2 /3 2 + |a cff | 2 ' 



(16) 

|beff| 2 /? 2 -|aeff| 2 
IM 2 /? 2 + kff| 2 ' 

Eq. ()15j) with R^ u replaced by computed as above but with a e fr(t) and 
6 e ff(t) replacing a e fr(t) and b e ff(t), respectively, gives the spin weight for 
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the events from B° decays. CP violating effects in the B° — > r + r~ and 
B° — > r + r _ decays are absent if Rq z = R z o, R xy = R yx and/or R xx = R xx . 

To simulate how the CP asymmetry in B^(B^) — > t + t~ decays are 
reflected in realistic observables the TAUOLA r-lepton decay library has been 
used (the method and numerical algorithm is given in ^2j). The input to 
the TAUDLA Universal Interface |13j is the spin density matrix of the 
t + t~ system resulting from the decay of a neutral particle. 

For simulations of the time integrated measurements, the time averaged 
matrix (R^ u ) has to be used 

m ,_ SdtV{Bl hYS {t)^r+r-)R, v {t) 
{ » V) - fdtT(B° hys (t)^T+r-) 

and {R^u) given by a similar formula. The asymmetry (JJJ) is then given by 

A i (i - (fl z0 ))r int - (l + (R z0 ))T int 

CP (l-(R g0 ))T int + (l + (R g0 ))r int ' 1 ' 

where 

r int =Jdt T(B° phys (t) - r+T"), f int = Jdt T(B° phys (t) - r+T"). 

AQp(ti, t-z) defined in (jSJ) is given by (fTTj) reversing the signs in the brackets. 

In the limit a = ±b{3 the two integrated rates r; nt and T- m i are equal 
resulting in (R z o) = — (Rzo)- It implies that even for 5c~p ^ there can 
be no CP violating effects in the observables sensitive to the longitudinal 
polarization of r's, like A^p, nor in R T . However, in this limit the ele- 
ments xx and xy of these matrices need not satisfy {R xy ) = — (R xy ) and 
(Rxx) = (Rxx)- Hence, the observables sensitive to transverse r polariza- 
tion can reveal CP violation even if the observables sensitive to longitudinal 
r polarization cannot. 



3. Supersymmetry scenarios 

In models of new physics Im(Ai) ^ Iu^A^ 1 ) and/or \q/p\ ^ 1 can be 
expected. However, if the decay rates are not at the same time strongly en- 
hanced, the detection of the fi + \i~ decay mode will become possible only at 
the LHC, while for the B® —> t + t~ mode at BELLE and BaBar the number 
of reconstructed events might be too small to detect any CP violation. 

Much more promising situation can occur in the supersymmetric sce- 
nario with a large ratio of the vacuum expectation values of the two Higgs 
doublets, v u /vd = tan (3 ~ 40 4- 50. Contributions from the Higgs penguin 
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diagrams with s-channel H° and A Higgs boson exchanges (31 IU d can 
become dominant easily saturating the experimental limits in Tabled if the 
Higgs particles H° and A in the range of order < 500 GeV even for other 
supersymmetric particle masses quite large, say in the TeV range [H I17j. 

In Ref. two different supersymmetric scenarios have been consid- 
ered: minimal (MFV) and non-minimal (NMFV) flavour violating, both 
with large ratio of VEVs. In both scenarios in which the B° ds -> l+l~ am- 




Fig. 1. The ratios R^ and R T as functions ofb/a for the phase Sqp — — ^arg(AL) 
0.1 (solid line), 0.3 (dashed), 0.5 (dotted) and 0.75 (dash- dotted). Ri(— Sqp) = 

plitudes have been dominated by the exchange of H° and A Higgs bosons, 
it was found that 

a ~ b or a ~ —b , 
(up to ^ 15%). For a = b the factors A^ and Xr simplify to: 
, Q 0*1- P q a* I + /3 

AL = ——5 , AR = " -z (18) 

p a 1 + p pal — p 

and all CP-sensitive quantities could be expressed in terms of one effective 
phase which can be taken as 

Scp = -^arg(A L ) (19) 

The effective CP phase is a function of the CKM phase and complex soft 
SUSY breaking Lagrangian parameters. 

The immediate consequence of a = b with \q/p\ ~ 1 is that for the \x~ 
final state, for which (3 = (1 — 4m^/M^) 1 / 2 is almost 1, the parameters \Xl\ 
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and \Xr\ assume values ~ 4 x 10 -4 and ~ 2.5 x 10 3 , respectively. As a result, 
the expected asymmetries are very small, at most \A^ P \ < 2 x 10~ 4 , |^4 CP | ~ 
10~ 3 . On the other hand, for the r + r~ final states /3 = (1 — 4m 2 /M^) 1 / 2 
differs substantially from 1 giving |A^| ~ 0.15, [A_r[ ~ 6.7 and the maximal 
possible values of the asymmetries are 

(A 1 CP ) max \=9% and |(A CP ) max | = 35% . (20) 

This is reflected in Fig. ^ in which possible CP violating effects in the ratio 
(jlljl for /z + n~ and t + t _ decay modes are shown as functions of b/a for four 
different values of the phase <5cp (keeping arg(a) =arg(6) and \p/q\ = 1). 
The plots demonstrate that the ratios R[ approach unity for a ~ ±£>/3, 
the feature which can be read also from the formula (jllj) if one takes into 
account that for a — > ±6/3 Ar ~ l/(a =F b[3) whereas |Ar| 2 ~ \a =F fy3| 2 - 
Therefore, for a ~ ±6 the deviation from unity of is tiny while for R T it 
can be quite substantial. 

The coefficients a and b (a and b) are constrained by the experimental 
limit in Tabled which in the case a ~ ±6 gives 

|a| ~ |6| < 4.9 x 10^ 9 . (21) 

In our simulations in the next section we conservatively set \a\ = \b\ < 
10~ 9 and treat both scenarios simultaneously, as all what matters are the 
values of \a\ = \b\ and the single CP violating phase 5cp which can be of 
order 1. As supersymmetric box and Z° penguin contributions as well as a 
finite difference of A and H° masses spoil the exact equality a = ±6, we 
investigate the effect of a and b different from each other by some 15 -i- 20%. 



4. Numerical results 

Two observables, known to provide valuable and complementary infor- 
mation on the spin state of decaying r lepton pairs, have been simulated: 

(a) 7r ± energy spectra in the decay channels r + — ► -k + v t (or t~ — > tt~v t ) 
reflect the longitudinal polarization of the individual leptons. Therefore 
they are sensitive to R z q and Ro z , as can be inferred from the expression 
(|15|) . i.e. to Re(a e fj6* ff ) as follows from (|16|). The CP violation is signaled 
if the energy spectrum of ir~(ir + ) originating from B°(B°) is different from 
the energy spectrum of 7r + (n~) originating from B (B ). In the following 
plots, the energy spectra are measured in the rest frame of the B°(B°) me- 
son assuming that the reconstruction of the event kinematics in the BELLE 
and BaBar experiments is sufficiently good for that purpose. 

(b) acoplanarity angle <p* between two planes spanned by the momenta 
of decay products of — > vr^vr coming from decays of both r leptons 
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into pv r ^Hl- It is sensitive to correlations between transverse components 
of r-lepton spins (i.e. to the elements R xx and R xy which in turn probe 
Im(a e ffb* ff ), as can be seen from (|16)0. The acoplanarity angle is defined as 

<P* = ( o ^ c l° T Sg i Pn ~ ' U+ \ < n (22) 
^ \ 27r - ^ for sgn(p„.- • n + ) > v ' 

where cos£ = an d two vectors n± = p^-i x p^o are normal to the 

planes determined by the momenta of pions which originate from p^ decays. 
Note that the full range of the variable < (p* < 2n of physical interest, 
and in addition events have to be sorted depending whether y±y2 > or 
2/12/2 < 0, where 

2/1 = p — 7TW~ ' m = w — TW~ ' ( 23 ^ 



since otherwise the spin correlations are washed out, as explained in Ref. [T| 
The acoplanarity distribution is evaluated in the rest frame of the p + p~ pair, 
but with the energies of 7r ± and 7r 0, s in (|23j) taken in the rest frame of the 
B°(B°). The difference between the distributions of the acoplanarity angle 
(p* measured in B° decays and the angle 2ir — cp* measured in B° decays for 
the same signs of yiy2 signals the CP violation. 

Fig. |21 shows the pion energy spectra and the acoplanarity distributions 
assuming \q/p\ = 1, a = b = and the CP violating phase 5cp = 

0.7. For all plots the same number of 5 x 10 5 r + r~ events from B^ and 
B^ decays has been generated with TAUOLA, although for the parameters 
chosen the ratio R T = 1.32, see Fig. ^ In the upper left panel the energy 
spectra of ir~ from B° decays (thick line) and of ir + from B° (thin line) are 
shown, while in the lower left panel shown are the spectra of 7r + from B 
decays (thick line) and of 7r~ from B° (thin line). The harder tt~ energy 
spectrum from B® decays than 7r + from B® in the upper left panel indicates 
that Br(B® — > t^Tr) > Br(B® — > t£t£), which is a clear signal of CP 
violation. In the acoplanarity plots (right panels) thick lines correspond 
to the distribution of (p* measured in B° decays, and the thin lines to the 
distribution of 2tt — if* measured in B° decays; in the upper right panel 
2/12/2 > 0, and y%y2 < in the lower right one. Different shapes of thick 
and thin lines seen in the right panels of Fig. [21 again indicate CP violation. 
In both energy and acoplanarity plots the CP violation is clearly seen and 
should be measurable even for small statistics. Note also that if upper and 
lower plots are combined, (i.e. no sorting according to the pion charge or 
sign of T/iT/2 is made), all CP asymmetries are lost. 

With decreasing [<5cp| the signal of CP violation deteriorates (especially 
in the pion spectra) and the possibility of distinguishing pion spectra and 
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0.25 0.50 0.75 1.00 1.0 2.0 3.0 4.0 5.0 6.0 



E„± /M B ip* 

Fig. 2. Results for the CP violating phase Sep = 0.7 and a = b. Left panels: Single 
7T^~ energy spectra in B°(B°) — > t + t~, t — > tt i/ T (i/ T ). In £/ie upper (lower) 
panel the thick line corresponds to the energy spectrum of tt~ ( of ir + ) from B° 
decays and the thin line to the energy spectrum of ir + (of tt~ ) from B° . Right 
panels: acoplanarity distributions of the p + p~ decay products in B°(B°) — > r + r _ . 
r — > j o ± ^ T (P r ), p —> n^ir . The thick lines correspond to the acoplanarity angle 
ip* measured in B° decays and the thin ones are for the angle 2ir — ip* measured in 
B° decays. Events in the upper (lower) panel have 2/12/2 > (2/12/2 <0). 

acoplanarity distributions from B° and B°, and hence the CP violation, 
would require increasingly large statistics which may not be attainable at 
BELLE and BaBar without major upgrades. 

As we discussed, the relation a = ±6 is only approximate. For b = 0.8 a 
with the same value of 5qp the CP violating effects in ir^ energy spectra get 
enhanced, while the acoplanarities are only slightly affected. On the other 
hand, for a approaching b(3 (for B — > t + t~ decays w 0.74) the effects of 
CP violation in the 7r ± energy spectra disappear, as expected and seen in 
the left panel of Fig. |31 In contrast, the acoplanarities shown in the right 
panel of Fig. |31 clearly indicate the CP violation even for a ~ b(3 confirming 
our earlier discussion demonstrating the complementarity of the energy and 
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acoplanarity distributions as a means to detect CP violation. 




E„± /M B ip' 

Fig. 3. As in the upper panels of figure^ but for Sqp — 0.7 and a = 0.8 b. 



5. Conclusions 

In the supersymmetry scenario with tan (3 ~ 40-^50 the rates of B®(B®) — > 
t + t~ decays are enhanced and could be detectable in the SLAC and KEK 
£?-factories. Moreover, the effective CP violating phase needs not be small. 
Therefore the CP asymmetries can be quite large as opposed to the B°(B°) — > 
decays in which they are kinematically suppressed. 

By using Monte-Carlo simulations we have investigated the possible ef- 
fects of CP violation in two realistic experimental observables and demon- 
strated that they might be detectable if the CP violating phase is reasonably 
large, i.e. 0(1). We have developed the necessary formalism and numerical 
tools allowing to apply the TAUOLA r-lepton decay library together with its 
universal interface to simulate fully the effects of the polarization of r + 
and t" originating from such decays. The tools can also be applied to de- 
termine the upper limits on the branching fraction of the B°(B°) — > r + r _ 
decays by the BaBar and BELLE collaborations. 
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